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The exact probability distribution of saturating states in random sequential adsorption 
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We consider the non-overlapping irreversible random sequential adsorption (RSA) process on 
one-dimensional finite line, which is known also as the car parking process. The probability of 
each coverage in saturating states is analytically and exactly obtained. In the derivation, a new 
representation of states in RSA process is introduced, which effectively works to make the calculation 
clear and simple. 
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I. INTRODUCTION 

We often see in various scientific and sociological fields 
that identical objects are sequentially placed at random 
on a region [l[ . Then, a meaningful question naturally 
arise: How much of the region is finally covered by them? 
That kind of phenomenon is known as a non-overlapping 
irreversible random sequential adsorption (RSA) , and the 
theories have been mainly developed in the field (see [l[ 
and references therein). It has been found that in ar- 
bitrary dimension, RSA process reach a saturating or 
jamming configuration, where further adsorption events 
are not possible. The final coverage as well as the tem- 
poral approach to the saturating state are of interest 
[H, i, [3,11, S, i, 0, H- Exact analytical results have been 
obtained mainly in one dimension, where the problem is 
also known as the parking problem [l], 0, S 3 • 

In those previous works, while the mean coverage has 
been obtained, the probability for each coverage in sat- 
urating states has never been revealed. Deriving it on 
one-dimensional RSA process analytically and exactly is 
the main purpose of this paper. 



II. THE MAIN RESULT 

Suppose identical particles of unit length arc placed 
sequentially at random on an interval of length I, subject 
to the constraint of no overlap. Here "place particles at 
random" means in this paper that each of them is placed 
based on the uniform distribution only on those empty 
spaces that have enough length to be placed. Note that 
it is not necessary to take the probability of finding 
already occupied place into account because we are not 
interested in the kinetics but only in the saturating 
states here. Then, the primary aim of this study is to 
obtain the probability of saturating states where there 
are n particles. The following is the main result, and its 
derivation is accomplished through this paper. 

In RSA process as above, the probability of saturat- 
ing states where there are n particles on the interval of 



length I becomes 
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where n C r '■= r \(n~r)\ an< ^ L"'J denotes the floor 
function. 



III. THE PROOF 

Let us introduce some notations which are used 
throughout this paper. Suppose now there are n par- 
ticles on the interval of length I. Then the state can be 
represented by a point in R" +1 . Namely, if we arrange 
the each length of n + 1 empty spaces between particles 
and between particles at the ends and the edges of the 
interval from the left side to the other, all of the states 
can be represented by a set of n+ 1 positive real number, 
say 

(oojffli, • • • ,a„); < a i} i = 0, 1, . . . ,n. (3) 
An example of this notation is shown in Fig. 1. 
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Fig. 1. The representation of states: The state 
illustrated above where n = 3 is expressed as 
(00,01,02,03) G K 4 . 
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Because the length of the interval is I, Oj (i 



0, 1, . . . , n) satisfies 



3=0 



(4) 



Thanks to the constraint condition ([5]), the number of 
independent valuables is reduced. Henceforth in this pa- 
per we eliminate clq. Thus, any state where there are n 
particles can be represented by a point in S n _i- n , 



(dl, . . ., O n ). 

Here Sn,p C R™ (0 < p) is defined as 



(5) 
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In S nj i- n , let us find the subset where all points cor- 
respond to saturating states. Those are such states 
that any length of n + 1 empty spaces are less than 
1. Therefore, the subset corresponds to R, l: i- n . Here 
Rn, P Cl" (0 < p) is defined as 



(ax,..., On) 



0<a t < 1 (i = l,... 
P - 1 < E"=i Oj < P 



■(7) 



On the other hand, all points in S nj i- n \R n ^- n satisfy 
1 < a, for some « = 0, 1, . . . , n. That is to say, states 
corresponding to them are not saturating states, and an- 
other particle can be placed on the interval. 

Next we discuss the probability of saturating states 
where there are n particles on the interval of length I. 
The derivation of the main result consists of three lem- 
mas. At first, we see the following lemma. 
Lemma 1. Suppose the particles are placed at random 
on the interval. Then, in terms of the notation (J3|), each 
state represented by a point in S n ,i-n has the same the 
probability. 

Proof. The proof is accomplished by mathematical in- 
duction. When n — 1, the event is placing a particle on 
the interval which is vacant. Therefore, the claim is ob- 
viously true by means of the definition of "at random" , 
which is definitely described before the statement of the 
main result. Assume the lemma is true when n = k 
for k = 1,2.... Then we can show the lemma is also 
true when n = k + 1 as below. Any state represented 
by a point in S k +i,i-(k+i) , sa Y (ai, ■ ■ ■ ,a k+ i) according 
to the notation ([5]), occurs if and only if a particle is 
placed in one of the k + 1 different states represented by 

(o 2 , 03, . . . , Ofe+l), (Oi + G 2 + 1,03, •••iGjfe+l)) (ai,02 + 

a 3 + 1, . . . , Ofe+i), (oi, 02, ■ • • , a k + a k+ i + 1), all 

of which are points in Sk,i—k- F° r example, a state 
(3/2, 2/5, 5/4) occurs if and only if a particle is placed in 
one of (2/5,5/4), (19/10,5/4) or (3/2,33/20). In short, 
any point in S^+i^-tk+i) is mapped from k + 1 different 
points in Sk.i-k by placing a particle. From the assump- 
tion of the lemma, these k + 1 mapping have the same 
probability. Therefore, it has been proved that each state 



represented by a point in Sk+i,i-(k+i) nas the same prob- 
ability if so does each state in S k ,i-k- Thus, it has been 
inductively proved the lemma holds for all n. □ 

Now we can calculate the probability of saturating 
states where there are n particles. In what follows, V(A) 
denotes the volume of domain A. Suppose there are n 
particles on the interval. Then, from Lemma 1, the con- 
ditional probability that this state is a saturating state 
in this given situation becomes 



V(RnJ-n) 
V(S n ,l-n) 



(8) 



On the other hand, the conditional probability that an- 
other particle can be added to this state becomes 



V(Sn,l-n\Rn,l-n) 
V(S n ,l-n) 



(9) 



Therefore, the probability of saturating states where 
there are n particles on the interval of length I, P/. n , 
becomes, 
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(10) 



What should be done to complete the derivation of the 



main result is to prove 



is equal to Qi. r in Eq. 



V(S„,,_ r ) 

(121). Then, let us calculate V(S n>p ) and V(R n , P ). These 
are shown in Lemma 2 and Lemma 3. 
Lemma 2. 



Proof. The domain S. 



n ,p 



P 

reads 



V{S n , p ) = ^r. 



S n ,p — ^ [Oi\ , ■ ■ ■ , fl n ) 

Therefore, the volume of Si iP becomes 
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= P -. □ (13) 
n\ 

By virtue of Lemma 2, let us find the volume of R n .p- 
Lemma 3. 



V{Rn, p ) = 



0. 



Z)iio(~i)' c 'i+i c 'i 
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(n<p - I). 
, (p-l<ri). 



(14) 



Proof. First of all, let us obtain the volume of 

. ,ri) 



T n , P := \ (a 



,a n ) 



< Oi < 1 (i = 1, . 

E"=l a 3 < P 



.(15) 
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Once we obtain V(T„ )P ), V(R n . p ) is calculated from 
V(Rn,p) = V(T n>p ) - V(T niP _i). To obtain the volume 
of T n>pi we subtract the superfluous domain from S ntP . 
In this way, it follows 



which corresponds to S n ,p-k- Therefore, since the Jaco- 
bian of the transformation is one, we obtain with the aid 
of Lemma 2 



V(T n , p ) = V(S n , p )- J2 V ( U n, 

l<ii <n 

+ --- + (-l) fe 

l<ii<...<i fc <n 
l<ii <...<i|_ p j <n 



E ^TC 2 ) 



ii— *LpJ 
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Here we have introduced for k = 1, . . . , |_pj , 



1 < a t {i = h, 
< a t (i ^ ii, 
E*=i a j < P 




Note that the formula (|T6|) is suitable only for [pj < 
that is p < n + 1. This is understandable if we concretely 
consider for n=3 as a simple example. For n + 1 < p, 
from the definition of T n>p (fToj) . it obviously follows 



1. 



(18) 



Let us calculate the right hand side of Eq. (fl6|) . For any 
combination of (ii, . . . , 



v(KT- lk ) = v(u% p -- k ) 



(19) 



because of the symmetry with respect to the permutation 
of upper indices. Then, it follows 

E ViKT Ak ) = nCuV{U% p k ). (20) 



Therefore, Eq. (fT6|) reads 



y(T n>p ) = Ec- 1 )*^^-*) (2i) 



,fe=i 



Let us obtain the volume of U^ p k . With the transforma- 



tion 



a, - 1, (i = 1, ...,k), 
at, (i = k + 1, . . . , n), 



(22) 



Un,p" k is written by 



(ai, . . . , a n ) 



< = ii, . . . , n) 

E"=i «i < p - fc 



(23) 



(p - k) n 



(24) 
(25) 



By substituting Eq. |(T3J) and Eq. flSSJ) into Eq. (|2Tj) . 
and according to Eq. (fl"8)) , we obtain 



(n<p-l), 
(p - 1 < ri). 



(26) 



Hence, we can immediately show that Eq. (|T4jl holds 
according to F(iJ„ iP ) = V(T ntP ) — V^T ntP -i), noting that 
V(T n>p ) = 1 for [p\ = n while Eq. Q15J) holds. □ 

As noted before, what must be done to complete the 



is 



derivation of the main result is to show that T 
equal to Qi^ r in Eq. ((2]). By substituting Eq. (fTTj) and 
Eq. (fblj) . We can immediately show that 



VCRn,p) / ' 
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(p-K 



IV. CONCLUDING REMARKS 

The RSA model on one-dimensional finite interval has 
been analyzed exactly. As a result, the probability of 
each coverage at saturating states has been obtained ex- 
plicitly for a generic system. Then we can obtain the ex- 
plicit expression of all moments. The main result, that is 
exact probability for any system size, would be useful for 
explaining the various actual RSA phenomena since ex- 
periments are always carried out on a system of finite size. 
Compared with the previous studies 0, 0, 0, 0, 0, 0, 0, Hi , 
this study takes another way to represent a state during 
the adsorbing process. Through the derivation of the 
main result, we see that this notation effectively works 
to make the calculation clear and simple. Besides this, 
omitting the kinetics also seems to be an important point 
to obtain the result. However, since we are almost always 
interested in not only final states but also kinetics in 
physics, construction of an appropriate master equation 
utilizing this result seems to be an significant problem 
for the future. The consistency with the previous studies 
[E 0, 0, 0, 0, 0, 0, Hi is another question. This should 
be investigated, if possible, by obtaining the asymptotic 
behavior of the mean coverage, that is the average of n/l, 
for I — > oo. 
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